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. ■ ABSTRACT 

^ , 

Oh, Context. The abundances of interstellar CH + and SH + are not well understood as their most likely formation channels are highly 
endothermic. Several mechanisms have been proposed to overcome the high activation barriers, including shocks, turbulence, and H 2 
vibrational excitation. 

Aims. Using data from the Herschel Space Observatory, we studied the formation of ions, in particular CH + and SH + in a typical high 
23 , UV-illumination warm and dense photon-dominated region (PDR), the Orion Bar. 

Methods. The HIFI instrument on board Herschel provides velocity-resolved line profiles of CH + 1-0 and 2-1 and three hyperfine 
transitions of SH + I2 - 0i. The PACS instrument provides information on the excitation and spatial distribution of CH + by extending 
' the observed CH + transitions up to / = 6 - 5. We compared the observed line intensities to the predictions of radiative transfer and 
PDR codes. 

qq , Results. All CH + , SH + , and CF + lines analyzed in this paper are seen in emission. The widths of the CH + 2-1 and 1-0 transitions 
■ are of ~5 km s , significantly broader than the typical width of dense gas tracers in the Orion Bar (~2-3 km s ) and are comparable 
' to the width of species that trace the interclump medium such as C + and HE The detected SH + transitions are narrower compared to 
( CH + and have line widths of ~3 km s -1 , indicating that SH + emission mainly originates in denser condensations. Non-LTE radiative 
transfer models show that electron collisions affect the excitation of CH + and SH + and that reactive collisions need to be taken into 
account to calculate the excitation of CH + . Comparison to PDR models shows that CH + and SH + are tracers of the warm surface 
region (A v <1.5) of the PDR with temperatures between 500 and 1000 K. We have also detected the 5-4 transition of CF + at a width 
of ~1.9 km s ', consistent with the width of dense gas tracers. The intensity of the CF + 5-4 transition is consistent with previous 
observations of lower- J transitions toward the Orion Bar. 

Conclusions. An analytic approximation and a numerical comparison to PDR models indicate that the internal vibrational energy of 
H 2 can explain the formation of CH + for typical physical conditions in the Orion Bar near the ionization front. The formation of SH + 
is also likely to be explained by H 2 vibrational excitation. The abundance ratios of CH + and SH + trace the destruction paths of these 
, ions, and indirectly, the ratios of H, H 2 , and electron abundances as a function of depth into the cloud. 



Key words. ISM: molecules - ISM: individual objects: Orion Bar 



1. Introduction 

The methylidyne cation CH + was one of the first molecules to 
be detected in the interstellar medium (Douglas & Herzberg 
1941). Early studies of CH + found its abundance to be consis- 
tently larger than the predictions of steady-state chemical mod- 
els in quiescent molecular clouds (e.g. Van Dishoeck & Black 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



1986). One of the possible formation routes is the endother- 
mic reaction C + + H2 + 0.41 eV — > CH + + H. To reproduce the 
observed CH + abundances, several mechanisms have been pro- 
posed to overcome the high activation barrier of the formation 
reaction. For low-density diffuse interstellar clouds, C-shocks 
(Pineau des Forets et al. 1986) and turbulent dissipation (Godard 
et al. 2009) have been proposed and confirmed by Falgarone et 
al. (2010a,b) and Godard et al. (2012). Alternatively, in denser 
regions with strong far-ultraviolet (FUV) radiation fields, the in- 
ternal energy available in the vibrationally excited H2 molecules 
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has been proposed to help overcome the large activation barrier 
(Sternberg & Dalgarno 1995, Agundez et al. 2010). 

Sulfanylium (SH + ) has a similar chemistry to CH + , hav- 
ing a formation route via S + and H2; however, this reaction is 
twice as endothermic as the CH + formation reaction. After non- 
detections of SH + in the UV domain in the spectra of nearby 
stars (Millar & Hobbs 1988, Magnani & Salzer 1989, 1991), 
the 526 GHz Nj = I2 - 0i transition of SH + has been detected 
in emission using Herschel toward the high-mass protostar W3 
IRS 5 (Benz et al. 2010). The 526 GHz transition has also 
been detected in absorption in the diffuse interstellar medium to- 
wards various distant star-forming regions (Godard et al. 2012). 
The 683 GHz transition of SH + has been detected in absorption 
towards Sgr B2(M) from the ground with the Carbon Hetero- 
dyne Array of the MPIfR (CHAMP+) receiver of the Atacama 
Pathfinder Experiment 12 m telescope (APEX) (Menten et al. 
2011). 

In this paper, we study the formation and excitation of CH + 
and SH + in a prototypical photon-dominated region (PDR), the 
Orion Bar. The Orion Bar is located at a distance of 414 pc 
(Menten et al. 2007). Its stratified structure has been the subject 
of many previous studies (such as Van der Wiel et al. 2009 and 
references therein). The mean density of the Orion Bar is about 
>10 5 cm 4 , the mean molecular gas temperature 85 K (Hoger- 
heijde et al. 1995), and the impinging radiation field is (1 - 4) x 
10 4 in Draine units (Draine field: y = 2.7 x 10~ 3 erg s _1 cm" 2 
for the energy range 6 < hv < 13.6 eV; Draine 1978). Most of 
the low-/ molecular line emission originates in an interclump 
medium with a density between a few 10 4 and 2xl0 5 ctrT 3 (Si- 
mon et al. 1997). High-density tracers such as HCN and H 13 CN 
originate in dense clumps, as confirmed by interferometric ob- 
servations (Young Owl et al. 2000). The density of the clumps 
is in the range between 1.5xl0 6 and 6xl0 6 crrT 3 (Lis & Schilke 
2003). Apart from the large clumps detected in H 13 CN deep 
inside the Bar, small, warm (7\in ~ 160 - 220 K), and dense 
(«h ~ 10 6 ~ 7 cm 4 ) condensations have been suggested to ex- 
plain the excited OH emission at the PDR surface (Goicoechea 
et al. 2011). 

This paper aims to characterize the medium where ions such 
as CH + and SH + form and to distinguish between the mecha- 
nisms that can overcome the high activation barriers of the for- 
mation reaction in a warm and dense PDR. We also report the 
detection of the CF + 5-4 transition. We will address another ion, 
OH + in a separate paper (Van der Tak & Nagy et al., in prepara- 
tion). 

2. Observations and Data reduction 

The CO + peak (a J20 oo = 05 h 35 m 20.6 s , feooo = -05°25'14") in 
the Orion Bar (Storzer et al. 1995) has been observed as part of 
the Herschel observations of EXtra-Ordinary Sources (HEXOS) 
guaranteed-time key program (Bergin et al. 2010) for the Het- 
erodyne Instrument for the Far-Infrared (HIFI, De Graauw et al. 
2010) of the Herschel Space Observatory (Pilbratt et al. 2010) 
in every HIFI band as a spectral scan. The data were reduced 
using the Herschel Interactive Processing Environment (HIPE, 
Ott et al. 2010) pipeline version 6.0. The velocity calibration of 
HIFI data is accurate to ~0.5 kms~' or better. The sideband de- 
convolution was done using the doDeconvolution task in HIPE. 
In this paper we use the HIFI bands la, lb, 2a and 3a from the 
HEXOS spectral line survey. These observations were carried 
out in March and April 2011 in load chop mode with a redun- 
dancy of 4 and had total integration times of 2.4 h, 2.2 h, 3.1 
h, and 1.3 h, respectively. The Wide-Band Spectrometer (WBS) 



backend was used which covers 4 GHz bandwidth in four 1 140 
MHz subbands at 1 . 1 MHz resolution. 

In addition to HIFI spectral scans observed as a part of the 
HEXOS key program, the CH + 2-1 transition was observed as a 
deep integration in a spectral scan in band 6b, with a total inte- 
gration time of 1 1 .7 h and a redundancy 4 in dual-beam-switch 
(DBS) mode. Both WBS and the High Resolution Spectrometer 
(HRS) backends were used. It was reduced using HIPE pipeline 
version 8.0. 

Besides the HIFI data, we used observations of the CH + 3-2, 
4-3, 5-4, and 6-5 transitions from the Photodetector Array Cam- 
era and Spectrometer (PACS, Poglitsch et al. 2010) onboard Her- 
schel. The PACS observations were carried out in September 
2010 and consist of two spectral scans in Range Spectroscopy 
mode with 5 range repetitions each (Joblin et al. 2012, in prepa- 
ration). The PACS spectrometer provides 25 spectra over a 
47" x 47" field-of-view resolved in 5x5 spatial pixels ("spax- 
els"), each with a size of ~9.4" in the sky. The measured width 
of the spectrometer point spread function (PSF) is relatively con- 
stant at A < 1 00 ^m, but it significantly increases above the spaxel 
size for longer wavelengths. The resolving power varies between 
A/AA -1000 (Rl grating order) and -5000 (B3A grating order). 
The central spaxel was centered at the same HIFI survey posi- 
tion. Observations were carried out in the "chop-nodded" mode 
with the largest chopper throw of 6 arcmin. The total integra- 
tion time was 3.2 h for the 1342204117 observation (B2B and 
Rl) and 2.7 h for the 13422041 18 observation (B3A). PACS data 
were processed using HIPE 6.0.3. 



CH + 3-2 co J = 6 - 5 - t r* (k) 




Fig. 1. PACS CH + 3-2 lines (black, centered on the CO + peak) overlaid 
on the distribution of the CO 7 = 6-5 peak brightness temperature 
(color image) observed with the CSO telescope at ~11" resolution (Lis 
et al. 1998). The PACS line intensity distributions are shown in units 
of Jy/spaxel as a function of wavelength in yum and are not velocity 
resolved. White contours show the brightest regions of v = 1-0 
S (1) emission (Walmsley et al. 2000). Lower-intensity W 2 extended 
emission is present in the entire field (Van der Werf et al. 1996). Violet 
stars show the position of the H I3 CN 7=1-0 clumps deeper inside 
the Bar (Lis & Schilke 2003). 



Table 1 shows the spectroscopic and observational parame- 
ters of the transitions used in this paper. The rest frequencies 
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are based on the Cologne Database for Molecular Spectroscopy 
(CDMS database, Miiller et al. 2005). In particular, the frequen- 
cies of the CH + 1-0 and 2-1 and 13 CH + 1-0 transitions are based 
on Miiller (2010). 

Table 2 includes the observed line parameters for our line 
sample observed with HIFI, which have been corrected for main 
beam efficiencies based on Roelfsema et al. (2012). In the case 
of the detected HIFI transitions we use the average of H- and 
V-polarizations. The detected lines show similar line profiles in 
both polarizations. The intensity difference between the polar- 
izations is 12% for CH + 7 = 1 - 0, 8% for CH + 7 = 2-1, 
-20% for SH + Nj = 1 2 - 0i, and -30% for CF + 5-4. To com- 
pare the line intensities of the transitions detected with HIFI with 
different beam sizes, we convert all the observed line intensities 
to a common ~47" resolution. We derive conversion factors be- 
tween the original beam sizes and ~47" based on the integrated 
intensity map of the HCN 4-3 transition from the James Clerk 
Maxwell Telescope (JCMT) Spectral Legacy Survey (Van der 
Wiel et al. 2009). For this, we assume the intensity distribution 
of the HCN 4-3 transition to trace the spatial structure of the 
Orion Bar. For the PACS data (Table 2), because of uncertain- 
ties in the PSF-correction, we use the mean value of the intensity 
measured at the central spaxel (/central) and that corresponding to 
the value integrated over 3x3 spaxels around the center (/3 X 3)- 
The corresponding error bars are the difference between / cen tiai 
and ^3x3. This method provides the correct agreement between 
Spectral and Photometric Imaging Receiver (SPIRE, Griffin et 
al. 2010) and HIFI data for CO lines (Joblin et al. 2012, in prepa- 
ration). The large (40%) error bar corresponding to the intensity 
of the 5-4 transition is caused by the presence of oscillations 
(fringes) in this range. Owing to a lower band intensity and a 
significant noise level, only /3 X 3 was measured for the intensity 
of the 6-5 transition; / ce ntrai was then estimated using a standard 
error value. 



3. Results 

3.1. The detected CH + , SH+, and CP transitions 

Figure 2 shows the velocity-resolved CH + transitions that were 
detected with significantly broad lines (~5 km s" 1 ) compared to 
the line width of dense gas tracers in the Orion Bar (~2-3 km 
s -1 , Hogerheijde et al. 1995 and an example of the CO 16-15 
transition in Fig. 2). The line width of the CH + 2-1 and 1-0 
transitions is not only significantly larger than the width of dense 
gas tracers, but also of species that trace a similar region to CH + , 
such as CO + (Sttirzer et al. 1995, Fuente et al. 2003). The 5 
km s _1 is comparable to the width of tracers of the interclump 
medium, such as C + (Ay ~ 3.8 km s , also shown in Fig. 2) and 
HF (Av ~ 4.9 km s 1 , Van der Tak et al. 2012). 

Other non velocity-resolved transitions from PACS, such as 
the 7=3-2 transition (Fig. 1), show extended emission detected 
over all PACS spaxels, decreasing with distance from the ion- 
ization front. The maximum CH + 3-2 emission is seen farther 
into the nearly edge-on PDR compared to the peak H2 u=l-0 
emission (Figure 1). As CH + forms via a reaction between C + 
and W 2 , we compare the CH + 3-2 intensity distribution to C + 
emission from Figure 4 in Ossenkopf et al. (2012). The C + peak 
matches the CH + peak within the HIFI beam. This indicates that 
CH + formation is limited by the C + abundance rather than by 
the H2 excitation that is traced by the H2 v = 1 — S( 1) emission 
observed by Walmsley et al. (2000). 

Figure 3 shows a tentative detection of the 13 CH + 7=1-0 
in V-polarization. If it were a real detection, the observed 



12 CH + / 13 CH + line ratio of ~40 would indicate an optical depth 
of the 12 CH + line of ~unity, assuming that the 13 CH + emission 
is optically thin. However, this seems unlikely since a 12 CH + 
optical depth of ~1 at a temperature >100 K (see Sect. 3.2) is 
inconsistent with the observed 12 CH + line intensity. Deeper ob- 
servations are needed to confirm the detection of 13 CH + . 

Figure 4 shows three hyperfine components of the Nj = I2 - 
0i transition of SH + (F = 3/2 - 1/2, F = 5/2 - 3/2, and F = 
3/2 - 3/2). The SH + F=3/2-3/2 transition is only detected in V- 
polarization. The line width of the detected SH + transitions (3.0 
km s -1 ) is narrower than those of CH + and are consistent with 
the width of dense gas tracers in the Orion Bar, suggesting that 
it does not originate in the same gas component as CH + . 

We have also detected the CF + 5-4 transition for the first 
time, with a line width of ~2 km s~'(Fig. 5). The 3-2, 2-1, 
and 1-0 transitions of CF + were previously detected toward the 
Orion Bar from the ground by Neufeld et al. (2006) with beam 
sizes of HPBW=24", 12", and 21", respectively. The velocity 
and the width of the 5-4 transition is consistent with the param- 
eters reported for the other detected transitions by Neufeld et al. 
(2006). One of the positions covered by Neufeld et al. (2006), 
05 h 35 m 22.8 s , -5°25'01", is close (ARA ~ 30", ADec- 13") 
to our observed position which is within the beam of HIFI at 
the frequency of the 5-4 transition (~44.2"). Assuming uniform 
beam-filling, a single excitation temperature of all four levels 
and optically thin lines, the measured line intensity is consistent 
with those detected by Neufeld et al. (2006) (Fig. 6) and implies 
a column density of ~2.1xl0 12 cirT 2 and a rotation temperature 
of -32 K. 

We derived upper limits for other non-detected CF + and SH + 
transitions. We estimated rms noise levels in the averaged spec- 
trum of H- and V-polarizations; 3<x upper limits on the integrated 
line intensities are estimated using (e.g. Coutens et al. 2012) 



7(3<x) [K km s" 1 ] = 3 rms V2 du FWHM, 

where dv is the channel width in km s and the rms noise level is 
derived in a velocity range of +5 km s -1 around the expected ve- 
locity. We use a full width at half maximum (FWHM) of 1.9 km 
s 'for the 13 CF + 5-4 (488664.3 MHz) and the CF + 6-5 (615365.6 
MHz) upper limits, and FWHM=3 km s^'for SH + Nj = 1 , - 0i 
(~683 GHz). The derived 3cr upper limits are listed in Table 2. 

3.2. Physical conditions traced by CH + and SH + 

To estimate molecular column densities, we use the non-LTE ra- 
diative transfer code RADEX (Van der Tak et al. 2007). We use 
H2 as a collision partner for the excitation of CH + , as we ex- 
pect a significant fraction of hydrogen to be in a molecular form 
at the observed position. Rates for inelastic collisions between 
CH + and H are not available, but are expected to be of the same 
order of magnitude as the rates for inelastic collisions between 
CH + and H2. If most H is in an atomic state, the density used 
as an input parameter is then the sum of n(H) and n(H2). We 
also include excitation via inelastic collisions between CH + and 
electrons, as the importance of excitation by electrons for HF 
has been recently demonstrated by Van der Tak et al. (2012). We 
apply an electron density of ~10 cm 3 . This is justified, if we 
assume that the electron abundance is determined by the abun- 
dance of C + . The column density of C + is approximately 10 18 
cirT 2 (Ossenkopf et al. 2012) and the H2 column density is ap- 
proximately 10 22 cirT 2 (e.g. Habart et al. 2010, Van der Wiel et 
al. 2009). This implies an electron abundance of 10~ 4 and using 
nQh) = 10 5 cirT 3 , an electron density of 10 cirT 3 . In the fol- 
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Table 1. Spectroscopic and observational parameters of the transitions used in this paper. 



Transition 










Frequency 

(MHZ ) 




A 

(s ) 




instrument/band 


beam-size 
( ) 


limb 


LH 1 — U 










83513 /.5 


A A 1 

40. 1 


6.36X10 


-3 


UTT7T U ~ « A 1 n 

Hlrl, band 3a 


26.5 


A 1 C 

0. /5 


CH + 2-1 










1669281.3 


120.2 


6.10x10 


-2 


HIFI, band 6b 


15.0 


0.72 


CH + 3-2 










2501440.5 


240.2 


2.20x10 


-1 


PACS 


9.4 1 




CH + 4-3 










3330629.7 


400.1 


5.38x10 


-1 


PACS 


9.4 1 




CH + 5-4 










4155872.0 


599.5 


1.07 




PACS 


9.4 1 




CH + 6-5 










4976201.4 


838.3 


1.86 




PACS 


9.4 1 




13 CH + 1-0 










830216.1 


39.9 


5.83x10 


.3 


HIFI, band 3a 


26.5 


0.75 


SH + Nj = 


h -0i 


, F = 


3/2- 


1/2 


526038.7 


25.3 


7.99x10 


-4 


HIFI, band la 


44.2 


0.76 


SH + Nj = 


h -0i 


, F = 


5/2- 


3/2 


526047.9 


25.3 


9.59x10 


_4 


HIFI, band la 


44.2 


0.76 


C I T 1 AT 

SH Nj = 


1 A 

I2 - Ui 


TP 


3/2- 


3/2 


jzolz4.y 


lj.5 


l.OUXlU 


-4 


Hlrl, band la 


A A O 

44. 1 


A 

U. /o 


CTT 1 \T _ 

on lyj — 


li - Ui 


TP — 
, F — 


J// - 


1 /o 


OojjjO. 1 


jZ.o 


Z.VUX 1U 


-4 


Hlrl, Dana za 


jj.Z 


U. / J 


SH + Nj = 


li-Oj 


, F = 


1/2- 


1/2 


683362.0 


32.8 


1.16x10" 


-3 


HIFI, band 2a 


33.2 


0.75 


SH + Nj = 


li-Oj 


, F = 


3/2- 


3/2 


683422.3 


32.8 


1.45x10" 


-3 


HIFI, band 2a 


33.2 


0.75 


SH + Nj = 


li-Oj 


, F = 


1/2- 


3/2 


683448.2 


32.8 


5.79x10 


-4 


HIFI, band 2a 


33.2 


0.75 


CF + 5-4 










512846.5 


73.8 


8.21x10 


-4 


HIFI, band la 


44.2 


0.76 


CF + 6-5 










615365.6 


103.4 


1.44x10 


-3 


HIFI, band lb 


44.2 


0.76 


13 CF + 5-4 










488664.3 


70.0 


7.10x10 


-4 


HIFI, band la 


44.2 


0.76 



1 The size of one PACS spaxel. 

Table 2. The detected CH + , SH + , and CF + transitions; the Gaussian fit parameters for the velocity-resolved transitions observed with HIFI and 
PACS; and upper limits for non-detections of other CF + and SH + transitions. The CH + J = 3 - 2, 4 - 3, 5 - 4, and 6-5 transitions are spectrally 
unresolved from PACS. The parameters from HIFI are based on the average spectrum of H- and V-polarizations, unless otherwise specified. 



Line 












J T MB dV 
(Kkms- 1 ) 


(km s- 1 ) 


AV 

(km s" 1 ) 


^peak 

(K) 


rms(r M B) 
(K) 


CH + 1-0 












24.9±0.20 


10.5+0.02 


5.46+0.04 


4.28+0.15 


0.12 


CH + 2-1 












10.6+0.20 


10.4+0.05 


4.57+0.11 


2.18±0.22 


0.23 


CH + 3-2 












2.19+0.31 










CH + 4-3 












1.01+0.19 










CH + 5-4 












0.42+0.15 










CH + 6-5 












0.17+0.06 










13 CH + 1-0 1 












0.46+0.11 


11.39+0.15 


1.30+0.34 


0.34+0.10 


0.15 


SH + Nj = 


1: 


-0! 


F 


= 3/2- 


1/2 


0.34+0.02 


10.92+0.09 


3.00 2 


0.11+0.02 


0.02 


SH + Nj = 


1: 


-0! 


F 


= 5/2- 


3/2 


0.57+0.02 


10.84+0.06 


3.00 2 


0.18+0.02 


0.02 


SH + Nj = 


1: 


-0! 


F 


= 3/2- 


3/2 


0.14+0.02 


10.34±0.24 


3.00 2 


0.04+0.01 


0.02 


SH + Nj = 


u 


-0! 


F 


= 3/2- 


1/2 


<0.18 






<0.15 3 




SH + Nj = 


h 


-0! 


F 


= 1/2- 


1/2 


<0.11 






<0.09 3 




SH + Nj = 


u 


-0! 


F 


= 3/2- 


3/2 


<0.13 






<0.12 3 




SH + Nj = 


h 


-0i 


F 


= 1/2- 


3/2 


<0.11 






<0.09 3 




CF + 5-4 












0.20+0.02 


11.13+0.11 


1.96+0.21 


0.10+0.01 


0.02 


CF + 6-5 












<0.08 






<0.08 3 




13 CF + 5-4 












<0.05 






<0.05 3 





1 Based on V-polarization data only. 

2 Fixed parameter in the fit. 

3 3xrms noise level. 



lowing we consider both H2 and electron collisions to probe the 
excitation of CH + and SH + . 

For CH + our calculations are based on collision rates from 
Turpin et al. (2010), for temperatures in the range between 
10 K and 200 K, covering transitions up to the 5-4 transition 
(F up = 599.5 K). These rates have been scaled from CH + -He to 
CH + -H2 based on Schoier et al. (2005). For electron collisions, 
we use collision rates from Lim et al. (1999) that are available 
for temperatures between 100 and 15000 K. 



Collisions of H2 and electrons are not always inelastic, but 
may lead to a chemical reaction. In the case of CH + this is im- 
portant since the collision rates with H2 and electrons are com- 
parable to the chemical reaction rates for CH + with H2 and elec- 
trons. For example, for CH + -e~ the chemical reaction rate is 
9xl0~ 8 cm 3 s~' for the destruction (Woodall et al. 2007) and is 
6.4xl0~ 7 cm 3 s~' for the excitation of the 7=1-0 transition 
at 1000 K. For CH + -H 2 the destruction rate is 1.2xl0~ 9 cm 3 s" 1 
(Woodall et al. 2007) and an excitation rate is 1 . 1 x 10~ 10 cm 3 s~' 
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Fig. 2. Top panel: Line profiles of CH + 2-1 and 1-0 transitions 
corresponding to the average of H- and V-polarizations observed with 
Herschel/HIFI toward the CO + peak in the Orion Bar. Bottom panel: 
Line profiles of C + and CO 16-15, for comparison, observed with Her- 
schel/HIFI toward the CO + peak in the Orion Bar. 
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Fig. 4. Line profiles of the three hyperfine transitions of SH + Nj = 
L - 0i observed with Herschel/HIFI corresponding to the average of H 
and V polarizations toward the CO + peak in the Orion Bar. 
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Fig. 5. Line profile of CF + 5-4 transition corresponding to the average 
Fig. 3. Line profile of 13 CH + 1-0 transition observed with Her- of H and V polarizations observed with Herschel/HIFI toward the CO + 



schel/HIFI in V-polarization toward the CO + peak in the Orion Bar. 



peak in the Orion Bar. 



for the / = 1 - transition at 100 K. Therefore, we consider the 
chemical formation and destruction rates in the statistical equi- 
librium calculation (e.g. Van der Tak et al. 2007). The statistical 
equilibrium for states i = 1 - N of energy that is solved using 
the RADEX code is given by the time-independent rate equa- 
tions 



dt 



-3 -1 

cm s , 



where 

Pij = Aij + BjjJ+Cij (Ej > Ej) 
= BijJ + dj (Ej < E j) 

and Ajj and Bjj are the Einstein coefficients, J is the mean in- 
tensity at the frequency of transition i — > j, C,j is the sum over 
all collision partners of the rates of inelastic, collision-induced 
transitions i — > j, rii is the number density (cm -3 ) of molecules 
in level i, and £), is the rate of destruction of the molecule in 
level i. When detailed knowledge of the state-specific formation 
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Fig. 6. Rotation diagram of CF + including the transitions detected by 
Neufeld et al. (2006). 



process is lacking, the formation rate into level i is expressed as 
a Boltzmann distribution over all states at an effective formation 
temperature Tj 

% oc giexpi-Ei/kTf), 

where gr, is the statistical weight of level i. When the destruction 
rate can be estimated, as for CH + here, then the total formation 
rate is normalized so that the total number density of molecules 
is consistent with its column density and the density of hydrogen 
in steady state. 

We assume that there is a balance between the formation 
and destruction of CH + . To simulate the chemical pumping 
effect described above, i.e. the effect of destruction and sub- 
sequent formation of CH + in excited levels, we add an arti- 
ficial level to the CH + level system, representing the dissoci- 
ated state, that is populated with a rate equivalent to the reac- 
tion rate of CH + -H 2 and CH+-e (Woodall et al. 2007). On 
the formation of CH + through the reaction of C + with vibra- 
tionally excited H2, the re-population from the dissociated level 
follows a Boltzmann distribution with a formation temperature 
of 7/ = 9920 K - 4560 K = 5360 K, where 4560 K is the 
required energy input for the endothermic CH + production and 
9920 K is the average energy of the vibrationally excited H2 lev- 
els following the 2-level approximation introduced by Rollig et 
al. (2006) where the full 15-level system was replaced by the en- 
ergetically equivalent 2-level system that provides the same total 
vibrational heating. 

Figure 7 shows the intensity predictions of two RADEX 
models with parameters in the range that can be expected for the 
Orion Bar. The error bars correspond to a 10% calibration error 
and a 10% error from obtaining integrated intensities by Gaus- 
sian fitting for the transitions observed with HIFI. The error bars 
corresponding to the PACS data are dominated by uncertainty 
in the PSF correction and are estimated as explained in Sect. 
2. Both electron collisions and H2 collisions suggest a kinetic 
temperature well above the average value (85 K) inferred for 
the interclump medium in the Orion Bar. Taking the formation 
pumping and collisional excitation described above into account, 
we find reasonable fits to the observed line-intensity distribution 
up to the 5-4 transition (the energy range for which the collision 
rates are available) with 7V(CH + )= 9 x 10 14 cirT 2 , T kin = 500 
K, and «(H2)=10 5 cm" 3 . The intensity of the 5-4 transition can 
be better reproduced with a kinetic temperature of 7\j n = 1000 
K. Temperatures between 500 K and 1000 K are expected near 



the edge of the cloud where the observations used in this paper 
have been taken. Assuming an electron density of n e = 10 cm" 3 , 
electron collisions mostly affect the two lowest-/ transitions for 
7\in = 500 K with an 10-13% increase in the intensities. 
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Fig. 7. Output of the RADEX models corresponding to a model with 
N(CH + )= 9 x 10 14 cm- 2 , n(H 2 )= 10 s cm 4 , n(e")= 10 cirT 3 , r kin = 
500 K. The blue symbols correspond to a model with excitation via 
collisions with H 2 , the green symbols to a model with excitation via H 2 
and electron collisions. 

In the case of SH + , calculations for collisions with H2 do not 
exist, so we use scaled radiative rates (J. Black, private com- 
munication) for a temperature range of 10-5000 K. We also 
use collision rates for electron-impact collisions calculated in 
the Coulomb-Born approximation (J. Black, private communi- 
cation), for a temperature range of 10-1000 K. Figure 8 shows 
the best fit models over-plotted on the observed line intensities. 
The error bars correspond to 20% of the observed line intensi- 
ties, including calibration error and the error introduced by the 
estimation of the integrated intensities using a Gaussian fit. We 
consider lower kinetic temperatures and higher volume densities 
than in the case of CH + , given that the line width suggests an 
origin from denser material. At the position of the CO + peak, 
warm (r kin ~ 160 - 220 K) and dense (10 6 " 7 cm" 3 ) conden- 
sations have been suggested to explain the OH emission. Us- 
ing SH + -H 2 collisions, a model with 10 6 cm" 3 , 7^ = 200 K, 
and 7V(SH + )=10 13 cm" 2 gives a reasonable fit to the observed 
line intensities (Fig. 8). On Fig. 8 we show a model with the 
same parameters, which includes electron collisions, assuming 
an electron density of « e = 10 cm" 3 . In both of these models the 
excitation temperatures are low (8.3-10.4 K) and the lines are 
optically thin (r ~ 0.02 - 0.2). 



4. The formation of CH + and SH + via H 2 vibrational 
excitation 

In this section we investigate the role of H2 vibrational excita- 
tion for the formation of CH + and SH + . We discuss alternative 
explanations in Sect. 6. 

4.1. Estimate based on an analytic approximation 

Testing if H2 vibrational excitation can drive CH + and SH + for- 
mation in an environment with a given radiation field and physi- 
cal parameters requires a detailed modeling with a PDR code, 
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Fig. 8. Output of the RADEX models for H 2 and electron colli- 
sions compared to the observed line intensities of SH + for a model 
with N(SH + ) = 10 13 cm" 2 , rc(H 2 ) = 10 6 cm 4 , n(e") = 10 cnr 3 , and 
T km = 200 K. 



with information on the chemical network and physical pro- 
cesses that affect the level populations of vibrationally excited 
Hh. A first indication can also be given using a simple analytic 
method to describe H2 vibrational heating with a two level ap- 
proximation (Rtillig et al. 2006). The H2 vibrational heating rate 
can be computed among all 15 vibrational levels in the ground 
electronic state, but neglecting the rotational structure, based on 
Equation C.2 in Rollig et al. (2006). A two-level system can be 
defined that results in the same vibrational heating rate as the 
full, 15-level system (Equation C.3, Rollig et al. 2006). The vi- 
brationally excited 'virtual' level has an upper level energy of 
A£eff=9920 K. The total rate for populating this vibrationally 
excited level is £0,1 = P\ X, where x is the radiation field in 
Draine units (Draine 1978) and P\ = 6.9 x 10" 10 s -1 is the 
formation rate of vibrationally excited H2 for the defined level 
of a radiation field of x = 1 ■ The de-excitation of this vibra- 
tionally excited level is via spontaneous emission, dissociation 
by the UV radiation field and by collisional de-excitation. The 
coefficient for spontaneous decay is A e ff = 1.9 X 10~ 6 s" 1 . The 
collisional de-excitation scales as n gas y e g with a rate coefficient 

of y e fr = 5.4 x 10~ 13 Vr s _1 crrT 3 . The dissociation rate is 
X x D e jf, where D e ff = 4.7 10 s -1 . These effective coefficients 
Aeff, Teff, and D e ff for the defined 2-level system as well as the 
energy of the defined vibrationally excited level A£eff are ob- 
tained by considering different asymptotic values of the density 
n and the radiation field x- By neglecting dissociation, the pop- 
ulation of the vibrationally excited level is dependent on the for- 
mation rate of vibrationally excited H2 as well as on the sponta- 
neous decay and collisional de-excitation rates. Calculating with 



igas = 10 s ctrT 3 and T = 500 K for the collisional de-excitation 
rate, there is a balance between these processes for a radiation 
field of x ~ 5 x 10 3 , which gives an expected lower limit on 
the radiation field, above which H2 vibrational excitation is ex- 
pected to be efficient enough to drive the formation of CH + and 
SH + . Even though neglecting dissociation introduces an addi- 
tional ~10% error, this calculation shows that for the radiation 
field in the Orion Bar (1 - 4 x 10 4 in Draine units), there is a 
large percentage of vibrationally excited H2 to react with C + and 
formCH + . This has been observed by Van der Werf et al. (1996) 
and Walmsley et al. (2000) and has already been noted for the 



formation of OH through the O + H2 — > OH + H reaction by 
Goicoechea et al. (2011). In the following section we test this 
idea with a more accurate approach, using PDR models. 

4.2. CH + formation 

We use the 1 .4.4 version of the Meudon PDR code (Le Petit et 
al. 2006, Goicoechea & Le Bourlot 2007, Le Bourlot et al. 2012) 
to model the observed CH + line intensities. This version in- 
cludes the Langmuir Hinshelwood and Eley-Rideal mechanisms 
to describe the formation of H2 on grain surfaces. The chemical 
pumping effect of destruction and formation on CH + level popu- 
lations is taken into account in addition to collisional excitation 
and de-excitation in the Meudon code (e.g. Gonzalez Garcia et 
al. 2008). The Meudon code treats CH + formation as described 
in Agundez et al. (2010) 



H 2 (j = - 7) + C + A CH + + H 



H 2 (v = 1) + 



CH + +H 



(1) 



(2) 



In reaction 1, H2 rotational levels up to J-l are used, which has 
an energy, £7=4586.4 K which is close to the activation barrier 
of the H2 + C + — » CH + + H reaction. We take into account the 
v — 1 vibrational level only because its energy (~5987 K) is 
enough to overcome the activation barrier of the CH + formation 
reaction. The formation rates are k\ = 1.58x 10~ 10 exp(-[4827- 
Ej/k]/T) based on Gerlich et al. (1987) and k 2 = 1.6 x 10~ 9 
(Hierl et al. 1997). 

We use isobaric models for typical conditions for the Orion 
Bar with pressures in the range between 5xl0 7 cirT 3 K and 
2xl0 8 cm 3 K, corresponding to T kin ~ 500 K (RADEX mod- 
els) and n~10 5 cm 3 ; and r kin ~ 1000 K (RADEX models) 
and n~2xl0 5 cirr 3 (typical interclump medium density, e.g. Si- 
mon et al. 1997); respectively. We apply a radiation field on 
the side where the cloud is illuminated from in the range be- 
tween ;ff r0 nt=l 4 and 3xl0 4 in Draine units (Draine 1978). We 
run the models up to a depth equivalent to a visual extinction 
of Ay ~ 10 mag. At the back side of the cloud (at Ay ~ 10 
mag), we use a radiation field 1000 times below that on the front 
Xback - Affront/ 1000. We adopt a cosmic-ray primary ionization 
rate of £ = 2 x 10~ 16 s per H2 molecule suitable for the dense 
ISM (Hollenbach et al. 2012). 

Figure 9 shows the results of a model for a pressure of 10 8 
cirr 3 K and a radiation field of^f ront =lxl0 4 , consistent with the 
radiation field near the ionization front of the Orion Bar. An in- 
clination of 60° was used to extract the line intensities, because 
of uncertainties in the computation of the line intensities in a ID 
model above this value (Gonzalez Garcia et al. 2008). This is 
reasonably close to the model with 75° inclination suggested to 
explain the geometry of the Bar (e.g. Melnick et al. 2012). The 
model with an inclination of 60° reproduces the observed CH + 
line intensities within a factor of 2 for the J=l-0 transition, and 
with an accuracy of 20% for the other transitions. Our RADEX 
calculations show the possible importance of electron collisions 
in the excitation of CH + . Therefore, to probe the effect of elec- 
tron collisions on the excitation of CH + , we implemented CH + - 
e~ collisions in the Meudon code. The models for P= 1 8 crrT 3 K 
and^ = 10 4 are shown in Fig. 9. Including electrons in the exci- 
tation of CH + , the model reproduces the observed line intensities 
with an accuracy of ~ 30%. Including electron collisions affects 
mostly the two lowest-/ transitions. The predicted intensity of 
the 7=1-0 transition increases by ~22%, and the intensity 
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of the 7 = 2-1 transition increases by ~18% after including 
electron collisions. 

The CH + abundance profile corresponding to this model is 
shown in Fig. 10 together with the gas temperature in the re- 
gion where CH + abundances peak. The model predicts that CH + 
forms near the surface of the PDR (Ay < 1) at high temperatures 
(T ~ 500 - 1000 K), consistent with the predictions by Agundez 
et al. (2010) and with our RADEX calculations. Though the best 
fitting models predict abundances to peak near the surface of the 
cloud at low Ay, the PACS observations of excited CH + used 
in this paper show a spatial extension along the area covered by 
PACS (47" x 47"), as shown in Fig. 1. SPIRE observations of 
the J = 1-0 transition (Naylor et al. 2010, Habart et al. 2010) 
show extended CH + emission over a ~ 200" x 200" region cen- 
tered on the crnooo = 05 h 35 m 22.83 s , S J2m = -05°24'57.67" 
position. The CH + J = 1 - emission mapped with HIFI was 
found to extend over a large region covering the OMC-1 cloud 
(Goicoechea et al., in preparation). One possibility is that the 
known dumpiness of the Orion Bar extends over a large volume 
and creates multiple PDR surfaces. Alternative explanation is 
that the extended CH + emission seen toward the region is the re- 
sult of a not completely edge-on PDR that is tilted to the line of 
sight. Models with lower pressures under-predict the observed 
line intensities. For example, a model with a pressure of 5xl0 7 
K cm" 3 underpredicts the line intensities with a factor of ~4. 
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Fig. 9. CH + line intensities as a function of E u , comparison between 
the observed line intensities and the predictions of an isobaric PDR 
model with a pressure of 10 s cm" 3 K for a radiation field of x = 10 4 - 



4.3. SH + formation 

SH + forms via a similar reaction to CH + , however, 
with an endothermicity about twice as high, AE=9860 K. 
Since state-to-state formation rates are not available for 
the H2 + S + — > SH + + H reaction, we use the rates of the 
H2 + C + — > CH + + H reaction as an approximation, taking into 
account H2 rotational levels up to £ 1 1 = 1 026 1 . 8 K and H2 in the 
v=l state up to £ = 10341.5 K. This is a reasonable assumption, 
since the total rates of the reactions for CH + and SH + formation 
via C + + H2 and S + + H2 are of the same order of magnitude 
(Woodall et al. 2007). To account for the higher activation bar- 
rier, we use k Umod = 1.58 x 10 10 exp(-[9860 - Ej/k]/T). We 
use the 1.4.4 version of the Meudon code (Le Petit et al. 2006, 
Goicoechea & Le Bourlot 2007, Le Bourlot et al. 2012), where 
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Fig. 10. Top panel: CH + and SH + abundances and the gas kinetic 
temperature as a function of A v for a pressure of P = 10 8 cm" 3 K and 
X = 10 4 . Bottom panel: H, H2, H2 (v=l), e", and C + densities. 

we introduce the SH + formation described above and use scaled 
radiative rates (J. Black, private communication) and electron- 
impact collisions calculated in the Coulomb-Bom approxima- 
tion (J. Black, private communication) for the excitation of SH + . 

With these assumptions, our best fit Meudon PDR model for 
CH + (P = 10 8 cm 3 K, x = 10 4 ) underpredicts the absolute 
intensities of the observed SH + transitions by a factor of ~3.5 
for the F = 5/2 — > 3/2 and F = 3/2 — > 1/2 transitions, and 
by a factor of 6.5 for the F = 3/2 — > 3/2 transition. Owing 
to the uncertainty in the formation rates, this agreement may be 
reasonable and suggests that like CH + , SH + can also be formed 
via H2 vibrational excitation in warm and dense PDRs. It may 
also suggest that SH + originates in a higher-pressure medium 
compared to CH + , which would explain the difference in the ob- 
served linewidths. The SH + abundances corresponding to this 
model are shown in Fig. 10. SH + abundances, like CH + abun- 
dances, peak near the surface of the cloud at Ay < 1 at high 
temperatures (500-1000 K). The SH+/CH+ abundance ratio in 
this region is between 0.01 and 0.1. This, however, is a lower 
limit on the SH + /CH + abundance ratio, since our model under- 
estimates the SH + line intensities. 



5. The destruction of CH + and SH + 

Figure 1 1 shows the abundance ratio of CH + and SH + predicted 
by our best-fit model in the region where CH + and SH + abun- 
dances peak. To understand these abundance ratios, it is essen- 
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Fig. 11. Abundance ratios of CH + and SH + in our best fit model 
(X = 10 4 , P = 10 s cnT 3 K) in the warm surface region (A v < 2) as a 
function of depth into the cloud. 

tial to study the destruction of CH + and SH + . CH + destruction 
can follow four main paths in the probed temperature and density 
regime 

H + CH + H 2 + C + 

, k H ,CH+ 

H 2 + CH + — — > H + CH+ 

e" + CH + -^-^ H + C 
/iv + CH + ^^C + H + 

Unlike CH + , SH + does not react with H 2 at the given physical 
conditions, as the reaction rate is orders of magnitude lower than 
that of reactions with electrons and H. Therefore the most impor- 
tant destruction paths are photodissociation as well as chemical 
reactions with H and electrons 

H + SH + H 2 + S + 

e~ + SH + — ^* H + S 
/ 1 y+SH + ^^H + S + 

The chemical reaction rates for these reactions are based on 
Woodall et al. (2007) and are summarized in Table 3. 



Table 3. Rates corresponding to the main destruction paths of CH + and 
SH + , based on Woodall et al. (2007). 



Reaction 


Rate 


Temperature regime 




(cm 3 s~') 


(K) 


&H.CH+ 


7.5 x 10~ lu 


10-41000 


£h 2 ,ch+ 


1.2 x 10~ 9 


10-41000 


k e -,CH + 


1.5 x 10~ 7 x (r/300)" ' 42 


10 - 300 


^H,SH+ 


1.10 x lO" 10 


10-41000 


£ e -,SH+ 


2.0 x 10~ 7 x (77300r ' 5 


10 - 300 


Reaction 


Rate 


Temperature regime 




(s" 1 ) 


(K) 


&hv,CH+ 


2.50 x 10- lu exp(-2.5A v ) 


10-41000 


^hv,SH+ 


3.00 x 10- 10 exp(-1.8A v ) 


10-41000 



At a depth equivalent to a visual extinction of Ay ~0.5, 
which is in the region where the CH + and SH + abundances 



peak, and at a temperature of ~830 K, most CH + (73.2%) is 
destroyed via collisions with H, while 24.3% is destroyed via 
collisions with H 2 and 2. 1 % is destroyed via collisions with elec- 
trons. Photodissociation is negligible in this regime, since it is 
responsible for 0.4% of the destructions of CH + . Deeper in the 
cloud, at a depth equivalent to Ay ~ 1, and at a gas temperature 
of ~570 K, most CH + (74.9%) is destroyed by collisions with 
H 2 , while small percentages of CH + are destroyed via collisions 
with H (22.8%), with electrons (2.2%), and by photodissociation 
(0.6%). 

At a depth of A v ~ 0.5, most SH + (-77.2%) is destroyed 
via reactions with H and a smaller percentage is destroyed 
via reactions with electrons (~18.2%) and by photodissociation 
(~4.4%). At a depth corresponding to Ay ~ 1, SH + is almost 
equally destroyed by reactions with H (~53.3%) and electrons 
(~43.7%). A smaller percentage of SH + is destroyed via pho- 
todissociation (~2.4%). 

At Ay ~ 1 and deeper, the SH + and CH + abundance ratio 
becomes higher than 0. 1 . Deeper in the cloud than Ay ~ 1 CH + 
abundances decrease more rapidly than SH + abundances. The 
[SH + ]/[CH + ] abundance ratio varies in the range between 0.01 
(at A v ~ 0.2) and >0.1 (at Ay ~ 1) and indicate that the de- 
struction of CH + becomes more efficient as a function of the 
depth into the cloud where H 2 takes over as the most important 
destruction partner, while with the decrease of atomic hydro- 
gen density SH + becomes more abundant because electrons are 
a less efficient destruction partner. While in diffuse clouds the 
abundance ratios of [SH + ]/[CH + ] trace the importance of shocks 
(e.g. Menten et al. 201 1) and properties of turbulent dissipation 
regions (e.g. Godard et al. 2012), in clouds exposed to high UV 
irradiation these abundances are sensitive to the abundance ra- 
tios of H, H 2 , and electrons as a function of depth into the cloud, 
which are determined by the radiation field strength of the irra- 
diation source. 



6. Discussion 

We have analyzed six rotational transitions of CH + and three 
transitions of SH + and reported the first detection of the / = 5 -4 
transition of CF + . We have shown that electron collisions affect 
the excitation of SH + and CH + , especially the lowest-/ transi- 
tions. We have also shown the importance of taking reactive col- 
lisions into account in the case of CH + excitation. We have con- 
firmed, both by an analytic approximation and by more detailed 
PDR modeling, that CH + formation is driven by H 2 vibrational 
excitation, unlike in the case of diffuse environments with lower 
UV radiation fields. SH + is also likely to form via H 2 vibra- 
tional excitation, although the lack of information on the exact 
state-to-state formation rates introduces an extra uncertainty in 
the models. 

6.1. The formation of CH + and SH + 

Spatially extended vibrationally excited H 2 emission was de- 
tected in the Orion Bar before the launch of Herschel (Van der 
Werf et al. (1996), Walmsley et al. (2000)), already indicating an 
importance in the chemistry of species that react with H 2 . Using 
Herschel, CH + was detected in the Orion Bar by Nay lor et al. 
(2010) and Habart et al. (2010), based on SPIRE maps of the 1-0 
transition. These observations show extended CH + 1-0 emission 
in the Orion Bar as well as in the OMC-1 cloud (Naylor et al. 
2010, Morris, P.; priv. comm., Goicoechea, J.; priv. comm.). 
Naylor et al. (2010) argue that the large spatial extent of CH + 
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into regions of low Ay suggests the importance of the forma- 
tion via H2 vibrational excitation. Our observations extend these 
studies, since the additional observed transitions up to 7=6-5 
provide additional evidence on the importance of the formation 
via vibrationally excited H2. An origin of CH + in the warm 
surface regions of the PDR is also confirmed by Goicoechea et 
al. (2011) who found a spatial correlation between excited OH 
2 n 3/2 / = 7/2" -> 5/2+ (-84.6 jum, observed with PACS) and 
CH + 3-2 emission, and that OH originates in the surface region 
(Ay < 1) of a high pressure gas component (10 s - 10 9 K cm -3 ). 
The formation and excitation of CH + in the Orion Bar is similar 
to that in the envelope of the high-mass protostar AFGL 2591, 
as its CH + emission can be explained to originate in the FUV- 
irradiated outflow-walls (Bruderer et al. 2010). Another region 
where CH + formation is driven by the strong FUV radiation field 
and can be explained by H2 vibrational excitation is the proto- 
planetary disc HD 100546 (Thi et al. 201 1), where CH + emission 
mostly originates in the outer disc and the disc surface in warm 
gas(r gas >400K). 

Other explanations for the formation of CH + applicable to 
the diffuse ISM include shocks (e.g. Pineau des Forets et al. 
1986). Tielens et al. (1993) have investigated that shocks do not 
contribute to the chemistry of the Orion Bar; therefore, we con- 
sider this scenario unlikely. Another scenario for CH + formation 
that has been successful in reproducing CH + abundances for the 
diffuse interstellar medium is the dissipation of turbulence (Go- 
dard et al. 2009, 2012). Though the CH + 1-0 and 2-1 transitions 
have broader line widths than most dense gas tracers in the Orion 
Bar, most of the lines detected in the Orion Bar are narrow (2- 
3 km s~'); therefore, we find it unlikely that turbulence plays a 
role in the chemistry of species detected in the Orion Bar. 

Unlike CH + , SH + has not been observed in a large variety 
of regions since its recent discovery in absorption toward Sagit- 
tarius B2 (Menten et al. 201 1). A recent study by Godard et al. 
(2012) probes CH + and SH + in absorption in the diffuse interstel- 
lar medium toward high-mass star-forming regions, suggesting 
a common origin for the formation and excitation of these ions, 
based on their observed linewidth-distributions and on compari- 
son with MHD shock models. However, in the diffuse ISM, SH + 
and CH + abundances are influenced by the dissipation of turbu- 
lence. SH + has also been detected in emission in the high-mass 
star-forming region W3 IRS5 (Benz et al. 2010), which repre- 
sents a region with physical conditions comparable to the Orion 
Bar, where the UV-radiation of the embedded protostars drives 
the chemistry of SH + . 

6.2. CH + and SH + as tracers of the warm PDR surface 

Though CH + and SH + most likely form via the same process and 
originate in the warm surface region of the PDR, a significant 
difference between CH + and SH + emission is suggested by the 
difference in the observed line widths. While the observed line 
width of SH + (At; ~ 3 km s~') is closer to that of dense gas 
tracers (Av ~ 2-3 km s~'), the width of the CH + 7=1-0 
and 2-1 transitions (Ai> ~ 5 km s~') is similar to that of HF 
(Av ~ 4.9 km s" 1 , Van der Tak et al. 2012) and C + (Au ~ 3.8 
km s _1 ), tracers of the interclump medium. The C91o' carbon 
recombination line was observed with the VLA with a width of 
2-2.5 km s (Wyrowski et al. 1997). It was found to match 
the H 2 [1-0 S(l)] distribution (Van der Werf et al. 1996) and 
its radial velocity was found to be consistent with that of H2 
pure rotational lines H 2 v=0-0 S(l), S(2), and S(4) (Allers et al. 
2005). The 13 C + lines have a slightly larger width of 2.5-2.8 km 
s , compared to that of the C91o' line. The larger width of the 



[Cn] 158 /im line compared to the C + lines can be a result of 
optical depth broadening of the C + line (with an optical depth of 
2-3). However, the C + line is also broader near the edge of the 
Bar, where the column density of material is lower, as are the 
line optical depths. In addition, the recombination line intensity 
is sensitive to the square of the electron density, while the fine 
structure line is sensitive to the local density only. Therefore the 
difference in line profiles of the C91o' and [Cn] 158 //m lines 
outside opacity broadening may also be related to gradients in 
the beam and along the line of sight, with the denser material 
having a lower velocity dispersion. 

This possible difference in the properties of the emitting re- 
gions is further indicated by our RAD EX models (Sect. 3.2). 
These models reproduce the observed CH + line intensities with a 
temperature of T = 500 - 1000 K and a density of n ~ 10 s cm -3 , 
but suggest a higher density component to explain SH + emission, 
T ~ 200 K, n ~ 10 6 cm" 3 , which is consistent with the properties 
of warm and dense condensations suggested to explain the ori- 
gin of excited OH (Goicoechea et al. 201 1) and high-/ CO line 
emission (Joblin et al. 2012, in preparation). In this case, ther- 
mal line broadening may contribute to the difference between 
the widths of the CH + and SH + lines. The expected contribution 

of thermal line broadening for CH + is Av—2 V21n2 J~^= 1 ■ 3 km 

s-'for 7\i n =500 K, and Au=1.8 km s 'for 7\ in =1000K. The con- 
tribution of thermal line broadening for SH + for 7\; n =200 K is 
Av-0.6 km s~ l . 

As an alternative explanation of the large observed line width 
of CH + , formation pumping may play a role in the broadening 
of CH + . As explained in Sect. 3.2, CH + formation results in an 
excess energy equivalent to 5360 K. This energy may be redis- 
tributed and go into kinetic motions. If the 5360 K excess energy 
goes into excess translational energy of the nascent CH + , and if 
this is identified as an ionic kinetic temperature upon formation, 
then the corresponding FWHM of Doppler motions is 4.4 km 
s- 1 . 

The difference between the widths of CH + and SH + may not 
originate in different excitation conditions, but in the difference 
in the chemistry of these ions. After its formation, CH + rapidly 
reacts with H and H2, therefore it is likely that its translational 
motions never become thermalized. In this case, the large ve- 
locity dispersion in CH + partially reflects the conditions of its 
formation. SH + on the other hand does not react rapidly with 
H (and H2), so that it is destroyed less rapidly by recombination 
with electrons. Therefore, SH + can become thermalized trans- 
lationally during its chemical lifetime. Therefore, while SH + 
traces the density and the temperature of the emitting region, 
CH + is more sensitive to the details of its formation process. 

Another interpretation of the broadening of CH + and other 
molecules could be that they originate in flows created by photo- 
evaporating clumps (e.g. Gorti & Hollenbach 2002, Mackey & 
Lim 2010). However, our observed data don't completely sup- 
port this assumption. If the FWHM of ions and other molecules 
detected at the same position had a contribution by the evaporat- 
ing flow, we would expect C + and CH + to have a similar flow ve- 
locity, as the momentum transfer from H2 in the H,+C + reaction 
is small. Based on our HIFI observations of the CO + peak, this 
does not apply, as FWHM(C + )~3.8 km s 1 and FWHM(CH + )~5 
km s . 

The difference in the width of reactive ions tracing the warm 
surface region of the PDR is a key part of understanding the 
chemistry of these ions. Future observations of the spatial dis- 
tribution of SH + will help to distinguish between these explana- 
tions. 



Article number, page 10 of 12 



Z. Nagy et al.: The chemistry of ions in the Orion Bar I. - CH + , SH + , and CF + 



6.3. An extension of the 'CF + ladder' in the Orion Bar 

Unlike CH + and SH + , CF + does not directly form via collisions 
with H2. The reaction between fluorine and H2 is followed by a 
reaction between C + and HF, where HF is the dominant reservoir 
of fluorine and was previously been detected in emission in the 
Orion Bar by Van der Tak et al. (2012) 

F + H 2 -> HF + H 

HF + C + -> CF + + H. 

CF + is the second most important fluorine reservoir and accounts 
for ~ 1 % of the gas-phase fluorine abundance. However, CF + has 
so far only been detected in two sources. The 2-1 and 1-0 rota- 
tional transitions have been recently detected with a spatially ex- 
tended emission toward the PDR in the Horsehead nebula (Guz- 
man et al. 2012). The first detection of CF + was toward the Orion 
Bar (Neufeld et al. 2006), showing spatially extended emission 
in the 1-0, 2-1, and 3-2 rotational transitions. Our observations 
extend the observed CF + transitions toward the Orion Bar up to 
the 5-4 transition and since the line intensity is consistent with 
the previously observed transitions, this work gives additional 
confirmation on the simple CF + chemistry tracing the surface 
layers exposed to UV irradiation. 

7. Conclusions and outlook 

We have analyzed six rotational transitions of CH + and three 
transitions of SH + and have reported the first detection of the 
5-4 transition of CF + . Our main conclusions are the following: 

- We have detected CH + up to the 6-5 transition. The 2-1 and 
1-0 transitions are spectrally resolved and show significantly 
broader lines (Av ~ 5 kms -1 ) than most dense gas tracers in 
the Orion Bar. SH + on the other hand shows significantly 
narrower lines than CH + (Av ~ 3 kms~'). Explanations of 
this difference include their origin in a different density (and 
temperature) component. Alternatively, because of its re- 
activity, CH + never becomes thermalized therefore, its ob- 
served properties trace the formation process rather than the 
properties of the emitting region, unlike for SH + . Informa- 
tion on the spatial distribution of SH + is needed to resolve 
this puzzle. 

- Inelastic collisions with H2 and electrons both affect the ex- 
citation of CH + and SH + , similar to the case of HF (Van der 
Tak et al. 2012). Reactive collisions are important in the ex- 
citation of CH + , but have less effect in the case of SH + . 

- Comparing the observed CH + intensities to predictions of 
PDR models for typical conditions in the Orion Bar, we 
confirm that CH + forms via reactions with vibrationally ex- 
cited H 2 , as predicted by Agundez et al. (2010). Our PDR 
models also show that CH + forms in the warm surface re- 
gion (T ~ 500 - 1000 K) of the PDR at high pressures 
(~10 8 cm' 3 K). 

- SH + is also likely to form via H2 vibrational excitation, as- 
suming that the formation rates are similar to that of CH + . 
SH + is also a tracer of the warm surface regions of the PDR. 

In the future, higher-resolution follow-up observations of a 
larger region in the Orion Bar will give more insight into the 
excitation conditions of SH + . Probing CH + and SH + formation 
in PDRs with a range of parameters, such as different radiation 
fields would help to deepen our understanding of the chemistry 
of these ions in regions exposed to UV irradiation. 
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